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PREFACE 


This  paper  is  directed  toward  potential  applications 
of  FNWC  spectral  ocean  wave  products  to  coastal  engineering. 
It  deals  specifically  with  the  formatting  and  properties 
of  climatological  wave  data  derived  from  synoptic  spectral 
wave  information,  and  speculates  or.  the  shoal  water  trans- 
formation of  these  deep-water  data. 

The  paper  was  presented  at  the  Fifteenth  Conference 
on  Coastal  Engineering  held  11-17  July  1976  in  Honolulu, 
Hawaii.  It  is  reprinted  here  with  permission  from  the 
1976  Conference  Proceedings,  Chapter  14  (pages  238-257) , 
published  1977  by  the  American  Society  of  Civil  Engineers. 
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CHAPTER  14 


OCEAN  WAVE  STATISTICS  FROM  FNWC  SPECTRAL  ANALYSES 

Warren  C.  Thompson,  Professor  of  Oceanography 
F.  Michael  Reynolds,  Lieutenant,  U.  S.  Navy 
Department  of  Oceanography,  Naval  Postgraduate  School 
Monterey,  California  93940 

ABSTRACT 

Climatological  wave  data  that  may  be  shoaled  and  refracted 
from  a deep-water  wave  station  can  be  compiled  in  two  forms  from 
spectral  ocean  wave  analyses  produced  by  the  Fleet  Numerical 
Weather  Central  at  Monterey,  California:  (1)  significant  wave 
statistics,  which  are  similar  to  statistical  tables  currently  in 
use,  and  (2)  spectral  element  statistics,  which  give  the  frequen- 
cy of  occurrence  of  energy  densities  contained  In  a matrix  of  15 
frequency  bands  and  12  direction  bands.  Experimental  formats  of 
both  types  of  statistical  compilations  are  presented,  their  pro- 
perties are  examined,  and  the  coastal  engineering  applications 
of  these  statistics  are  discussed. 


BACKGROUND 

For  certain  purposes  coastal  engineers  require  f requency-of - 
occurrence  statistics  on  ocean  waves  compiled  for  a selected  deep- 
water wave  station  located  offshore  from  the  coastal  sector  of 
interest.  These  wave  climatology  tables  characteristically  give 
the  statistical  frequency  with  which  hel ght- perl od-dl recti  on 
( H - T - ip ) data  elements  occur  per  average  month  or  average  year. 

They  are  of  particular  value  In  beach  erosion  studies  because  of 
the  fact  that  the  H - T - elements  can  be  shoaled  and  refracted  from 
deep  water  to  any  location  in  Intermediate  or  shallow  depths  up  to 
the  breaker  point.  The  gradient  of  the  statistical  wave  field  in 
deep  water  is  generally  sufficiently  small  along  the  coast  that 
wave  statistics  prepared  for  a given  deep-water  wave  station  may 
be  carried  into  shoal  water  depths  for  many  miles  In  either  direc- 
tion from  the  station. 

Frequency-of-occurrence  statistics  are  derived  from  synoptic 
weather  maps  by  the  wave-hi ndcast  method.  Heretofore,  the  analy- 
sis of  weather  maps  for  wave  generation  parameters  (wind  velocity, 
fetch,  and  wind  duration)  has  been  performed  manually.  Because 
of  the  large  amount  of  labor  Involved,  climatological  data  of  this 
type  are  not  available  for  most  coasts  of  the  open  oceans.  Sta- 
tistical tables  have  been  compiled  by  the  hindcast  method  for  the 
Atlantic,  Gulf,  and  Pacific  coasts  of  the  United  States,  however, 
and  these  have  had  wide  use  (tables  currently  In  use  were  prepared 
by  Saville,  1954;  Newmann  and  James,  1955;  Bretschnei der  and  Gaul, 
1956;  Marine  Advisers,  1961;  and  National  Marine  Consultants,  1960, 
1961).  The  effort  Involved  in  making  manual  hlndcasts  has  also 
limited  these  compilations  to  three  years  of  weather  maps.  The 
yearly  variability  of  the  statistical  wave  field  at  a given  open 
ocean  wave  station  In  middle  and  higher  latitudes  Is  known  to  be 
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appreciable;  accordingly,  a three-year  data  base  for  long-term 
average  statistics  must  be  considered  minimal. 

In  December  1974  the  Fleet  Numerical  Weather  Central  (FNWC) 
of  the  United  States  Navy  at  Monterey,  California  put  into  routine 
operation  a spectral  ocean  wave  computer  model  that  produces  "real- 
time" computations,  termed  analyses,  of  the  deep-water  wave  field 
in  spectral  form  12-hourly  for  a grid-point  field  covering  the 
major  oceans  of  the  Northern  Hemisphere  (waves  are  not  computed 
for  the  Southern  Hemisphere).  The  wave  field  is  computed  in  3- 
hour  time  steps  from  the  surface  wind  field,  which  in  turn  is 
derived  from  the  observed  sea-level  pressure  field.  The  grid-point 
spacing  is  approximately  150  nautical  miles.  Details  concerning 
the  FNWC  Spectral  Ocean  Wave  Model,  including  comparisons  between 
FNWC  wave  computations  and  open-ocean  observations,  are  given  by 
Lazanoff  and  Stevenson  (1975). 

In  addition  to  issuing  daily  wave  analyses  and  forecasts,  FNWC 
has  begun  running  historical  weather  maps  through  their  spectral 
wave  program  and  is  in  the  process  of  storing  the  6-hourly  wave 
field  for  the  Northern  Hemisphere  oceans  on  magnetic  tapes.  Pres- 
ent plans  call  for  the  production  of  20  years  of  6-hourly  analyses 
for  the  period  1955-1974.  When  the  project  is  completed,  which 
will  require  many  months,  the  archived  wave  data  will  provide  the 
potential  to  extend  the  present  wave  statistics  coverage  geograph- 
ically to  most  coasts  and  islands  of  the  open  North  Pacific,  North 
Atlantic,  and  North  Indian  Oceans,  and  to  Increase  the  time  base 
nearly  an  order  of  magnitude  as  well. 

Wave  statistics  for  coastal  engineering  use  have  not  previous- 
ly been  produced  from  FNWC  spectral  wave  computations.  In  this 
paper, forma ts  for  presenting  f requency-of-occurrence  statistics 
are  illustrated  and  applications  of  the  statistics  are  discussed. 
Additional  information  is  given  by  Reynolds  (1976). 


DEVELOPMENT  OF  THE  TABLES 

The  approach  taken  by  the  authors  to  develop  wave-occurrence 
statistics  from  FNWC  spectral  analyses  was  to  mock  up  suitable  for- 
mats manually  from  a short  series  of  FNWC  computer  printouts  for 
a given  wave  station.  This  permitted  experimentation  with  the  for- 
mat design,  including  the  selection  of  suitable  scales.  The  wave 
station  chosen  is  a FNWC  grid  point  located  at  latitude  50.9°N  and 
longitude  145. 6°W  near  Ocean  Weather  Station  Papa  In  the  open  Gulf 
of  Alaska  (a  deep-water  site  along  an  open  ocean  coast  could  have 
been  chosen).  A summer  month  and  a winter  month,  August  1974  and 
February  1975,  were  selected  for  extremes  in  seasonal  wave  activ- 
ity and  with  the  expectation  of  a dominance  of  swell  during  one 
month  and  sea  in  the  other.  The  statistical  tables  for  February 
only  are  presented  here.  The  equivalent  tables  for  August  can  be 
found  in  Reynolds  (1976).  It  may  be  noted  that  the  August  wave 
statistics  reveal  much  lower  wave  activity,  shorter  wave  periods. 
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and  a higher  incidence  of  swell  versus  sea  than  occurred  in  Feb- 
ruary. Because  of  missing  synoptic  analyses,  the  data  for  Febru- 
ary that  are  shown  in  the  tables  were  drawn  from  45  12-hourly 
computer  printouts  totalling  22*s  days  of  time. 


SYNOPTIC  SOURCE  DATA 

Familiarity  with  the  FNWC  wave  analysis  printout  is  required 
for  a full  understanding  of  the  design  and  characteristics  of  the 
tables  presented  herein.  A sample  printout  for  the  selected  wave 
station  for  00Z  on  31  August  1974  is  shown  in  Table  1. 

The  15  wave-frequency  values  shown  across  the  top  of  the  print- 
out (Item  1)  each  represent  the  central  frequency  of  a given  FNWC 
frequency  band.  The  central  frequencies  are  seen  to  cover  a period 
range  from  6.1  seconds  (0.164  Hz)  to  25.7  seconds  (0.039  Hz).  The 
limits  and  width  of  each  frequency  band  are  given  in  Table  2.  The 
12  wave-direction  values  shown  in  the  right-hand  column  in  degrees 
(Item  2)  and  in  the  left-hand  column  in  coded  form  (Item  3)  each 
represent  the  central  direction  of  a 30°  directional  band.  The 
geographical  directions  are  seen  to  be  irregular.  They  vary  from 
grid  point  to  grid  point  because  of  the  global  map  projection  used 
by  FNWC  (Lazanoff  and  Stevenson,  1975). 

The  block  of  180  numbers  forming  the  body  of  the  table  (Item  4) 
constitutes  the  f requency-di  rection  (f-ip)  spectrum  of  the  wave 
field.  Each  number  represents  the  wave  energy  density,  expressed 
as  variance  ( f t ^ ) , contained  in  the  frequency  and  direction  bands 
indicated.  The  sums  of  the  col umns  ( I tem  5)  represent  the  wave 
energy  contained  in  the  15  frequency  bands  irrespective  of  wave 
direction  and  constitute  the  familiar  f-spectrum.  The  sums  of 
the  rows  (Item  6)  give  the  wave  energy  contained  in  the  12  direc- 
tion bands  irrespective  of  wave  frequency,  and  collectively  repre- 
sent the  less  familiar  ^-spectrum. 

The  sums  of  the  component  variance  values  contained  in  the 
f-ip  spectrum,  the  f-spectrum,  and  the  ^-spectrum  (Items  4,  5,  and 
6)  are  identical,  and  for  the  waves  being  considered  are  seen  to 
total  3.145  ft’  (Item  7).  It  is  of  interest  to  note  that  two 
well-defined  wave  trains  from  different  sources  are  revealed  in 
the  printout.  The  larger  set  is  a swell  train  with  a dominant 
period  of  about  15.0  seconds  (0.067  Hz  band)  and  a dominant  direc- 
tion from  the  West  (259°  band);  the  other  set  is  sea  of  about  9.7 
second  period  (0.103  Hz)  from  the  South  (199°). 

At  this  point,  some  comment  is  needed  regarding  the  wave  fre- 
quency bands  (Item  1)  represented  in  the  printout.  The  15  bands 
have  been  carefully  selected  by  FNWC  to  minimize  computing  time 
and  yet  co'/er  the  full  frequency  range  for  ordinary  gravity  waves 
with  suitable  bandwidths.  It  may  be  noted  In  Table  2 that  the 
bands  below  a frequency  of  0.080  Hz  are  of  constant  width  (0.0055 
Hz),  but  that  at  higher  frequencies  the  bandwidth  has  been  widened 
by  multiples  of  2,  3,  4,  and  18.  A consequence  of  this  design  is 
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that  in  order  to  properly  compare  the  magnitude  of  two  energy 
density  values,  the  values  must  be  adjusted  to  a common  frequency 
bandwidth.  Choosing  the  common  or  standard  bandwidth  to  be  Af  = 
0.0055  Hz,  this  can  be  done  by  multiplying  the  energy  density 
values  in  the  printout  by  the  appropriate  energy  density  factors 
listed  in  Column  4 of  the  table.  By  way  of  example,  the  second- 
ary peak,  of  the  f-spectrum  (Item  5 of  Table  1)  appears  to  lie  at 
a frequency  of  about  0.133  Hz,  whereas  this  maximum  actually  lies 
at  about  0.103  Hz  (the  primary  peak  of  the  f-spectrum  lies  at 
about  0.067  Hz).  The  energy  content  of  these  two  frequency  bands 
are  seen  from  the  printout  to  have  values  of  0.41  ft2  per  Af  = 

0.0165  Hz  (from  Table  2)  and  0.33  ft2  per  Af  = 0.0110  Hz,  or 
about  0.14  ft2  and  0.16  ft?  respectively,  when  adjusted  to  the 
common  bandwidth  of  Af  = 0.0055  Hz.  The  component  variance  values 
of  the  f-i i>  spectrum  (Item  4)  must  be  similarly  adjusted  for  direct 
comparison  of  their  energy  contents. 

It  should  be  further  noted  that  If  the  variance  values  compos- 
ing the  f-spectrum  shown  in  the  printout  (Item  5)  were  plotted  on 
a linear  frequency  scale,  the  resulting  graph  would  be  exaggerated 
at  the  high-frequency  end  of  the  spectrum  because  of  the  choice  of 
variable  frequency  bandwidths.  Nevertheless,  the  sum  of  these 
variance  values,  3.145  ft2,  correctly  represents  the  total  energy 
content  computed  for  the  waves  present.  For  a correct  graphical 
representation  of  the  frequency  spectrum  (to  the  accuracy  permitted 
by  the  widths  of  the  frequency  bands),  the  variance  values  must 
first  be  adjusted  to  a common  frequency  bandwidth  through  use  of 
the  energy  density  factors.  No  such  difficulty  arises  with  the 
^-spectrum  because  of  the  use  of  equal  directional  bandwidths. 

The  FNWC  decision  to  represent  the  frequency  spectrum  by  unequal 
bandwidths  affects  the  way  in  which  statistical  tables  may  be 
formatted  and  the  statistics  compiled,  and  requires  care  by  the 
user  in  order  to  assure  correct  interpretation  of  the  data.  This 
is  commented  upon  further  in  the  discussion  of  the  individual 
tables. 

Returning  to  the  spectral  printout  in  Table  1,  the  waves  com- 
puted to  be  present  at  the  station  at  the  designated  time  are  seen 
to  be  described  in  considerable  detail  by  180  AV-Af-At elements, 
each  being  a component,  AV,  of  the  total  variance  of  the  wave 
field  (V)  contained  in  a specific  frequency  band,  Af,  and  direc- 
tion band,  A ip . These  wave  data  are  referred  to  herein  as 
spectral  e 1 ement  data . The  wave  information  in  the  printout  may 
also  be  summarized  by  designating  the  significant  height, 
the  period  of  maximum  energy  density,  T , and  the  ' 

direction  of  maximum  energy  density,  ? For  the, wave  field 
shown  in_the  printout,  fl,,.,  = 6.8  feetm  (R,,.,  = 4V^  where  V = 

3.145  ft^),  T lies  in  /Jthe  frequency  ' band  centered  at 
15.0  seconds  m x(or  0.67  Hz),  and  lies  in  the  direction  band 

centered  at  259  degrees.  Thus,  them  *wave  field  can  be  described 
in  a gross  way  by  a single  Ri / Q-Tmay"^mav  data  unit.  Wave  data 
in  this  form  will  be  ' referred  to  as 

significant  wave  data. 
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It  is  evident  that  wave  statistics  derived  from  a time  series 
of  FNWC  spectral  analyses  may  be  compiled  as  (1)  significant  wave 
stati  sti  cs , consisting  of  one  fli/o-T  -1(1  data  unit  per  FNWC 
analysis , and  as  (2)  spectral  max  ax  e 1 ement  statistics, 

composed  of  180  AV-Af-Aij/  data  units  per  FNWC  analysis  [or  fewer 
than  180  units  when  variance  values  of  0.00  ft*  are  excluded 
from  the  statistical  tabulation).  Example  formats  of  each  are 
given  in  the  following  sections. 

In  the  statistical  tables  and  the  discussion  herein,  wave  fre- 
quency and  wave  period  are  used  interchangeably,  in  a reciprocal 
sense,  according  to  convention;  frequency  is  used  when  reference 
is  made  to  the  FNWC  spectral  model  and  period  when  operational 
wave  products  are  of  interest. 


SIGNIFICANT  WAVE  TABLE 

The  significant  wave  statistics  compiled  for  the  selected 
station  for  February  1975  from  FNWC  spectral  wave  analyses  are 
presented  in  Table  4.  Because  of  the  sparcity  of  the  data  in 
the  individual  direction  bands,  the  entries  for  all  12  directions 
have  been  combined  into  a single  table  for  illustration  purposes 
in  this  paper;  tables  intended  for  practical  use  would  display 
each  direction  component  separately.  The  values  in  the  table  re- 
present the  number  of  12-hourly  synoptic  analyses  occurring 
during  the  221-s  day  data  period.  The  height  scale,  coded  for 
brevity  (code  given  in  Table  3),  is  linear  with  an  increment  of 
2 feet.  Regarding  the  period  scale,  wave  occurrences  were  com- 
piled from  the  printouts  by  the  FNWC  frequency  bands  shown  in 
Table  2,  but  for  purposes  of  convention  each  column  of  Table  4 
is  labelled  for  the  period  equivalent  to  the  central  frequency 
of  each  FNWC  frequency  band.  The  equivalent  period  bandwidths 
are  seen  from  Table  2 to  be  variable  and  to  lie  mostly  in  the 
range  of  1 to  2 seconds. 

The  design  of  Table  4,  except  for  the  use  of  frequency  as 
a base  instead  of  period,  is  essentially  identical  to  existing 
f requency-of-occurrence  tables,  but  with  one  Important  excep- 
tion. In  compiling  wave  statistics  derived  from  the  manual 
analysis  of  weather  maps,  wherein  wave  trains  generated  in  in- 
dividual wind  areas  are  computed  and  propagated  separately,  it 
has  been  the  general  practice  to  record  each  set  of  waves  as 
an  independent  occurrence  whenever  two  or  more  trains  are  pres- 
ent simultaneously  at  the  wave  station.  This  results  In  tables 
containing  more  hours  of  wave  occurrence  In  a month  than  actual- 
ly occurs,  particularly  when  separate  tables  for  sea  and  swell 
are  compiled.  This  presents  the  user  with  the  problem  of  how  to 
adjust  the  data  to  100%  occurrence.  In  addition,  higher  wave 
heights  must  occur  in  nature  than  these  tables  generally  indi- 
cate because  multiple  wave  trains  are  a common  occurrence  at  open 
ocean  wave  stations.  In  the  FNWC  Spectral  Ocean  Wave  Model,  the 
growth  of  sea  and  the  propagation  of  swell  are  complexly  inter- 
related and  individual  wave  trains  are  not  differentiated  in  the 
f-if/  spectrum.  Accordingly,  the  statistics  compiled  from  the 


FNWC  spectral  analyses  reflect  wave  heights  associated  with  the 
statistical  occurrence  of  multiple  wave  trains,  and  the  super- 
occurrence problem  does  not  exist. 

In  using  significant  wave  tables  compiled  from  FNWC  spectral 
analyses,  it  may  be  desirable  for  some  purposes  to  distinguish 
swell  from  sea  and  to  estimate  the  relative  age  of  the  swell, 
i.e.,  whether  it  has  travelled  over  a short,  moderate,  or  long 
decay  distance.  Because  swell  steepness  diminishes  with  increas- 
ing travel  distance  from  the  generating  area,  some  measure  of 
the  steepness  of  the  wave  field,  y , can  be  used  for  this  purpose. 
In  this  study,  y is  defined  in  terms  of  H,/3  and  T by  analogy 
to  monochromati c waves  as  shown  in  Table  5.'JThe  xrel ati onsh i p 
selected  between  wave  age  and  ys  is  given  in  the  lower  part  of 
the  table.  Thus,  young  swell  has  a statistical  steepness  in  the 
range  of  about  1/40  and  1/100,  has  travelled  0 to  250  nautical 
miles  from  the  location  where  it  was  generated,  and  has  been  re- 
duced in  significant  height  by  a factor  of  1.0  to  0.5.  It  may 
be  observed  that  because  significant  wave  tables  are  a matrix  of 
values  of  R.|  ,3  and  Tmax»  they  are  also  tables  of  y and  accord- 
ingly of  wave  age.  This  permits  wave-age  boundaries  to  be 

superimposed,  as  in  Table  4. 

It  may  be  recalled  that  the  term  "wave  age"  was  introduced 
by  drup  and  Munk  (1947)  to  describe  the  growth  stages  of  sea 
wave  generation.  They  defined  it  as  the  ratio  of  wave 
surface  wind  speed,  C/U,  and  related  it  to  wave  steep- 
• n extending  use  of  the  term  to  swell,  which  has  consider- 
practical  merit,  wave  age  is  defined  differently  than  for 
sea  but  remains  related  to  the  steepness  of  the  wave  field. 

The  zone  labelled  "sea"  in  Table  4 contains  all  synoptic  oc- 
currences where  sea  is  the  dominant  component  of  the  wave  field 
present,  regardless  of  the  wind  speed  generating  the  sea  or  of 
its  stage  of  growth.  The  zone  is  seen  to  have  a cutoff  such 
that  wave  heights  greater  than  a certain  value  do  not  occur  for 
a given  period.  Young  seas,  because  they  are  of  maximum  steep- 
ness, should  lie  along  this  data  envelope,  whereas  fully  arisen 
seas  should  fall  near  the  inner  margin  adjacent  to  the  zone 
labelled  "young  swell".  A well-defined  data  cutoff  should  be 
expected  whenever  a significant  number  of  synoptic  entries  repre- 
sent seas,  but  may  not  be  evident  in  statistics  where  swell  is 
dominant. 


SPECTRAL  ELEMENT  TABLES 

Spectral  element  statistics  compiled  for  the  selected  station 
for  February  1975  from  FNWC  synoptic  wave  analyses  are  shown  in 
Table  6.  The  entries  in  the  table  were  read  directly  from  the 
180  component  variance  values  contained  In  each  of  the  45  12- 
hourly  printouts  available  for  the  month  (eg..  Item  4 of  Table  1). 
The  entries  for  all  12  direction  bands  have  been  combined  into  a 
single  table  for  purposes  of  comparison  with  Table  4. 


Table  6 is  identical  to  Table  4 with  respect  to  the  period 
scale  and  the  manner  of  representing  frequency  of  occurrence 
(i.e.,  by  the  number  of  12-hourly  analyses);  however,  the  signi- 
ficant-height scale  is  replaced  by  a variance  scale  chosen  to 
give  a reasonable  spread  to  the  data.  The  variance  scale  is 
coded  (Table  3)  to  simplify  the  table.  should  be  noted  that 

the  units  of  variance  in  Table  6 are  ftVFNWC  frequency  band- 
width. 

Tables  6 and  4 were  derived  from  the  same  FNWC  printouts  and 
the  total  energy  content  represented  by  their  entries  is  identi- 
cal. Nevertheless,  there  are  no  values  common  to  both  distribu- 
tions, and  consequently,  there  are  no  means  by  which  the  data  in 
the  two  tables  can  be  equated  numerically.  Neither  can  one 
table  be  derived  from  the  other.  These  circumstances  arise  from 
the  fact  that  although  the  energy  in  the  waves  at  a given  place 
and  time  can  be  separated  analytically  into  spectral  components 
and  can  subsequently  be  reconstituted  to  yield  a significant  wave 
height  representing  that  sea  surface,  the  same  spectral  components 
when  distributed  in  a statistical  table  with  components  from 
other  synoptic  analyses  can  no  longer  be  identified  and  re- 
assembled to  produce  the  significant  height. 

A comparison  of  these  two  tables  offers  some  interesting  con- 
trasts, however,  and  gives  some  feeling  for  the  coarseness  of 
significant  wave  data  as  a descriptor  of  ocean  wave  climatology. 
The  spectral  element  table  contains  some  fifty  times  more 
data  than  the  significant  height  table;  specifically,  2126 
variance  entries  (>  0 f 1 2 ) versus  45  height  entries.  The  variance 
entries  are  distributed  in  a wider  range  of  period  bands,  from 
6.1  to  25.7  seconds  versus  9.7  to  18.0  seconds,  and  extend  into 
what  are  ordinarily  considered  very  long  periods.  Further,  the 
entry  of  maximum  energy  is  centered  in  the  20.0  second  band  for 
the  spectral  element  statistics  versus  the  16.4  second  band  for 
the  significant  wave  statistics.  It  should  additionally  be  noted 
that  the  wave  activity  in  February  1975  was  appreciable;  this  is 
indicated  in  Table  4 by  the  small  number  of  occurrences  of  waves 
under  10  feet  (Code  05  and  under),  but  is  not  apparent  in  Table  6 
where  most  of  the  entries  are  In  the  lowest  energy  categories. 

A peculiarity  of  Table  6 is  that  a given  energy  category  or 
code  number  does  not  represent  the  same  energy  density  for  all 
frequency  bands.  This  is  due  to  the  variable  widths  of  the  FNWC 
frequency  bands.  For  example,  Code  03  represents  an  energy  level 
of  0. 50-0. 74„ft2/o. 0055  Hz  bandwidth  in  the  12.0  second  band  but 
0.50-0.74  ftV0.0110  Hz  bandwidth  In  the  adjacent  10.9  second 
band,  giving  the  10.9  second  band  about  half  the  energy  density 
of  the  12.0  second  band.  Accordingly,  energy  density  values 
cannot  readily  be  compared  across  the  table.  This  circumstance 
wou  1 d compl i ca te  the  use  of  this  type  of  table  and  would  likely 
result  in  incorrect  interpretation  of  the  statistics. 

This  difficulty  may  be  avoided  by  reducing  the  tomponent 
variance  values  to  the  standard  frequency  bandwidth,  Af  * 0.0055 


Hz,  through  multiplication  by  the  energy  density  factors  given 
in  Table  2.  This  operation  may  be  applied  either  to  the  180 
variance  values  in  each  FNWC  printout,  or  with  a slight  reduction 
in  statistical  accuracy  to  the  statistics  compiled  directly  from 
these  values.  Table  7 was  derived  from  Table  6 by  the  latter  2 
method.  The  variance  scale  of  Table  7 thus  has  the  units  of  ft  / 
0.0055  Hz  standard  bandwidth  for  all  periods  in  the  table. 

Table  7 has  the  desirable  quality  that  a given  energy  code 
number  represents  the  same  energy  level  for  all  frequency  bands. 
On  the  other  hand,  the  table  contains  less  total  energy  than 
Table  6.  This  is  due  to  the  fact  that  in  adjusting  the  eneigy 
of  the  FNWC  frequency  bands  to  the  standard  bandwidth  in  the 
periods  from  6.1  to  10.9  seconds,  segments  of  the  frequency  spec- 
trum lying  outside  the  0.0055  Hz  bandwidth  centered  about  each 
period  are  in  effect  excluded. 

A feature  of  the  spectral  element  tables  of  particular  in- 
terest that  is  best  seen  in  Table  7 is  the  well-defined  energy 
cutoff  extending  diagonally  across  the  table.  The  cutoff  repre- 
sents the  energy-saturated  frequencies  of  sea  spectra  associated 
with  the  strongest  winds  occurring  during  the  month.  The  cause 
and  nature  of  the  cutoff  can  be  more  readily  visualized  from 
cons i dera t i on  of  Table  8 which  shows  the  spectral  entries  such 
as  might  be  extracted  from  a single  FNWC  printout  for  a fully 
arisen  sea  of  50  knots  (in  the  absence  of  swell).  The  entries 
in  this  table  show  the  energy  density  per  0.0055  Hz  and  30°  band- 
widths  (as  in  Table  7)  for  the  condition  in  which  the  mean  wind 
direction  coincides  with  the  center  of  the  direction  band  (this 
condition  gives  the  maximum  variance  values  possible  for  the 
50-knot  sea).  The  reader  may  recognize  this  distribution  as  a 
frequency  spectrum  graph  in  tabular  form  that  has  been  rotated 
180°  about  both  the  energy  density  and  frequency  axes.  If  one 
now  imagines  a graph  of  nested  frequency  spectra  for  fully  arisen 
seas  generated  over  a range  of  wind  speeds,  the  spectra  are  seen 
in  the  higher  frequencies  to  occupy  a relatively  narrow  band  on 
the  graph,  with  the  outermost  spectrum  being  associated  with 
the  strongest  wind  and  representing  the  energy  cutoff  in  the 
equivalent  statistical  tables.  If  non-fully  developed  spectra 
for  the  same  wind  speeds  are  added  to  the  imagined  graph,  the 
unsaturated  portions  of  each  will  lie  inside  this  band,  as  will 
also  the  saturated  low-frequency  tails  of  fully  arisen  spectra 
associated  with  weaker  winds. 

It  is  evident  from  these  considerations  that  the  entries  in 
Table  7 lying  inside  the  cutoff  zone  may  represent  (1)  saturated 
frequency  components  associated  with  fully  arisen  seas  produced 
by  weaker  winds,  (2)  unsaturated  components  of  seas  produced  by 
strong  or  weak  winds,  or  (3)  swell.  There  is  no  way  to  associate 
these  entries  with  sea  or  swell  however,  whether  or  not  this 
might  be  desired  for  practical  purposes.  All  that  can  be  said 
in  this  vein  is  that  the  entries  lying  In  a band  along  the  data 
envelope  clearly  represent  dominantly  sea.  It  should  be 
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mentioned  here  that  the  energy  cutoff  in  Table  7 is  not  fully 
analogous  to  the  height  cutoff  in  Table  4;  occurrence  entries 
rep  resenting  sea  may  be  found  anywhere  inside  the  data 
envelope  of  Table  7 but  lie  only  in  the  zone  labelled  "sea" 
i n Table  4. 

The  wave  statistics  for  February  are  presented  in  still  an- 
other form  in  Figure  1.  This  figure  was  constructed  by  cumulat- 
ing the  occurrences  in  Table  7 from  high  to  low  energy  density 
within  each  period  band  (cumulative  table  not  shown),  and 
plotting  the  cumulated  values  against  the  period.  The  curves  in 
the  graph  are  isopleths  of  energy  density  coded  according  to 
Table  3.  Figure  1 combines  the  advantage  of  Table  6 in  that  the 
total  energy  represented  is  correct,  and  Table  7 in  which  the 
energy  densities  are  referenced  to  the  standard  0.0055  hz  band- 
width. This  graphical  presentation  gives  the  user  a distinctly 
better  visual  impression  of  the  statistical  distribution  than 
do  the  two  tables;  however  tables  directly  provide  the  numerical 
values  a user  needs. 

Figure  1 is  of  particular  interest  not  so  much  as  a format 
in  which  spectral  element  statistics  may  be  presented  but  because 
of  the  possibility  which  it  suggests  for  describing  climatologi- 
cal wave  data  by  a mathematical  function,  at  least  for  „ome  ocean 
regions  and  seasons.  The  energy-occurrence  curves  in  vbe  figure 
are  seen  to  resemble  nested  wave  frequency  spectra  for  fully 
arisen  seas  associated  with  a range  of  wind  speeds,  with  the 
dashed  curve  being  the  counterpart  of  the  T curve  for  the 
spectra.  The  fact  that  the  pattern  of  thism  x small-sample  dis- 
tribution is  so  well  defined  suggests  that  the  long-term  average 
pattern  for  February  is  similar  at  this  station.  It  may  be 
expected,  on  the  basis  of  continuity,  that  similar  wave  climatology 
patterns  occur  in  adjacent  months  and  over  adjacent  regions  of 
the  North  Pacific  Ocean;  accordingly,  it  seems  possible  that  all 
of  these  distributions  might  be  described  by  a function  having 
variable  coefficients  that  provide  for  seasonal  and  geographical 
gradients.  From  a practical  standpoint,  the  availability  of 
wave  climatology  models  for  ocean  regions  would  greatly  improve 
data  packaging,  allow  greater  versatility  In  data  manipulation, 
and  improve  the  prediction  of  infrequently  occurring  wave  condi- 
tions. 

Consideration  of  the  spectral  characteristics  of  sea  and  swell 
further  suggests  that  in  ocean  regions  and  seasons  where  swell  is 
dominant,  the  wave  climatology  may  not  be  amenable  to  mathematical 
description.  A graph  similar  to  Figure  1 prepared  for  August  (not 
shown)  when  the  swell-to-sea  ratio  Is  substantially  higher  at  the 
selected  station  is  inconclusive  on  this  point  because  of  the  low 
frequency  of  wave  occurrence  in  that  month.  It  should  be  recalled 
that  Figure  1 is  a data  composite  of  all  12  direction  bands.  If 
long-term  average  statistics  for  February  were  displayed  for  each 
direction  Land,  it  is  reasonable  to  expect  that  a pattern  similar 
to  that  in  the  figure  would  appear  for  all  or  some  of  the  component 
di rections . 
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USES  OF  THE  STATISTICS 


Significant  wave  statistics  tables  derived  from  FNWC  spectral 
analyses  are  similar  in  construction  to  climatology  tables  current- 
ly in  use,  except  for  tabulation  of  the  wave  data  by  variable 
period  increments.  This  design  feature  may  appear  untidy  but 
should  not  restrict  the  uses  of  the  statistics. 

The  coastal  engineering  applications  which  would  be  made  of 
these  tables  appear  to  be  identical  to  those  of  existing  wave  cli- 
matology tables.  Reasons  for  favoring  use  of  the  FNWC-derived 
significant  wave  statistics,  however,  are  that  they  are  free  of 
the  familiar  super-occurrence  problem  in  which  tables  contain  more 
total  hours  of  wave  occurrence  than  actually  occurs,  and  that  they 
take  into  account  simultaneous  occurrence  of  wave  trains  arriving 
at  a wave  station  and  therefore  contain  higher  wave  heights. 

With  regard  to  spectral  element  statistics,  there  are  con- 
siderable uncertainties  concerning  their  uses  and  apparently  no 
experience  available  to  provide  guidance.  The  following  comments 
are  therefore  limited  to  the  more  obvious  things  that  can  and 
cannot  be  done  with  these  statistics. 

The  AV-Af-Aii)  values  composing  spectral  element  tables  compiled 
for  a deep-water  wave  station  can  be  shoaled  and  refracted  in  the 
same  manner  as  the  spectrum  of  a given  set  of  waves  (described  by 
Pierson,  Neumann,  and  James,  1955);  thus,  the  spectral  element  cli- 
matology can  be  determined  for  a selected  shoal-water  site  (except 
near  the  surf  zone  since  the  breaker  height  cannot  be  identified 
in  the  spectral  statistics).  These  statistics  cannot  be  converted 
into  wave  height  statistics,  however,  and  so  cannot  be  used  to 
compute  design  wave  heights,  wave  forces,  breaker-height  statistics, 
or  1 i ttora 1 -dr i f t rates  by  current  methodology;  these  quantities 
must  still  be  calculated  using  the  more  gross  significant  wave 
statistics. 

The  principal  uses  for  wave  statistics  in  spectral  element 
form,  both  in  coastal  waters  and  the  open  ocean,  appear  to  lie  in 
the  realm  of  resonant  reaction  of  structures  to  waves.  Practical 
applications  can  readily  be  visualized  involving  piling  platforms, 
spar  platforms,  tethered  breakwaters,  floating  causeways,  surface 
vessels,  and  surface-effects  craft.  The  nature  of  the  transform 
from  spectral  wave  information  to  structures  response  and  how  the 
transform  can  be  applied  quantitatively  appear  to  be  little  under- 
stood, however. 

Although  specific  uses  for  spectral  element  climatology  are 
not  now  developed,  the  marked  differences  which  these  statistics 
display  from  the  significant  wave  statistics  derived  from  the 
same  synoptic  source  data,  along  with  the  possibility  of  modelling 
these  spectral  distributions  mathematically,  compels  further 
attention  to  these  data. 
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Table  2:  FNWC  FREQUENCY/PERIOD  PARAMETERS 
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Table  3:  CODE  FOR  EXPRESSING  SIGNIFICANT  WAVE  HEIGHT  AND 

SPECTRAL  ENERGY  DENSITY  IN  THE  CLIMATOLOGY  TABLES 

Significant  wave  table  Spectral  element  tables 
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5.50 

- 

5.99 

17 

32.0 

- 

33.9 

17 

6.00 

- 

6.49 

18 

34.0 

- 

35.9 

18 

6.50 

- 

6.99 

19 

36.0 

- 

37.9 

19 

7.00 

- 

7.49 

20 

38.0 

- 

39.9 

20 

7.50 

- 

7.99 

21 

40.0 

- 

41.9 

21 

8.00 

- 

8.49 

22 

42.0 

- 

43.9 

22 

8.50 

- 

8.99 
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ALL  DIRECTIONS 


FEBRUARY  1975 


CENTRAL  PERIOD  (SECONDS) 

6.1  I 6 5 I 7.5  I 8 6 I 9 7 1 10.9  |12.0 1 12.9 1 13.8 |15.0 116.4 118.0 1 20.0 1 22. 5 125. 7 
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Table  5:  DEFINITIONS  OF  WAVE  STEEPNESS  AND  WAVE  AGE 
(Applicable  to  deep  water  only) 

WAVE  STEEPNESS,  y 


Monocbromati c 

waves 

rm 

H 

H H * 

wave  height 

5.12  T2  L * 

wave  length 

T « 

wave  period 

Spectrum  w.  • 

„ = hn 

fll/3  fll/3 

* significant 
wave  height 

* period  of  max 
energy  density 

A Tmax  5* 

: AGE  (in  terms  of  y$) 

12  T2 

max  T 

max 

Wave  age 

Ys 

Significant 

height 

reduction* 

(H0/HF) 

Decay 
di stance* 

( naut  mi ) 

Sea 

Young  swel  1 

1/40 

1.00 

0 

Moderate  swell 

1/100 

0.50 

250 

Old  swell 

1/250 

0.25 

1600 

HF  = sea  height  in  the  generating 

area 

Hp  = swell 

height 

at  decay  distance  D 

‘Approximate  values  for  swell  generated  in  an  average 
extratropi cal  cyclonic  storm  from  the  SMB  wave  fore- 
casting graphs  (Bretschneider,  1958). 
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ALL  DIRECTIONS 


FEBRUARY  1975 


Q 

UJ 

Q 

O 

O 


9 

5 

o 

z 

< 

ffi 

O 

LU 

a 


i 

z 

u. 

■>». 

> 

K 

(/> 

z 

UJ 

a 

> 

a 

cr 

UJ 

z 

Ul 


CENTRAL  PERIOD  (SECONDS) 


6.1 

6.5 

75 

86 

9 7 

10.9 

12.0 

129 

13.8 

15.0 

16  4 

180 

o 

o 

CM 

22  5 

25.7 

01 

142 

135 

122 

123 

123 

153 

143 

158 

129 

108 

76 

60 

37 

11 

13 

02 

141 

115 

137 

40 

41 

33 

39 

37 

29 

27 

17 

2 

1 

03 

1 

65 

49 

33 

30 

17 

17 

16 

9 

9 

8 

1 

i 

04 

65 

51 

18 

19 

16 

D 

D 

6 

3 

H 

1 

i 

05 

2 

29 

19 

20 

13 

13 

9 

6 

1 

1 

06 

20 

9 

8 

D 

a 

D 

1 

07 

18 

4 

8 

4 

2 

1 

■ 

08 

16 

3 

6 

3 

2 

■ 

09 

1 

3 

3 

a 

1 

1 

3 

10 

5 

■ 

2 

11 

6 

. 1 

3 

12 

■ 

5 

1 

1 1 

1 

13 

1 

2 

1 

14 

15 

1 

1 

16 

1 

17 

18 

19 

1 

20 

21 

1 

22 

Table  6:  SPECTRAL  ELEMENT  STATISTICS 

Energy  density  (coded)  per  FNWC  frequency  bandwidth 
AND  30  DIRECTION  BANDWIDTH 


ENERGY  DE 


ALL  DIRECTIONS 


FEBRUARY  1975 


CENTRAL  PERIOD  (SECONDS) 


01  283 


02 


10.9 

12.0 

12.9 

13.8 

15.0 

186 

143 

158 

129 

108 

- 1 

51 

39 

37 

29 

27 

39 

30 

17 

17 

16 

34 

19 

16 

HD 

20 

13 

13 

9 

76 

60 

37 

11 

17 

12 

2 

1 

9 

9 

8 

1 

6 

3 

1 

B 


1 1 


1 


Table  7:  SPECTRAL  ELEMENT  STATISTICS 
Energy  density  (coded)  per  0.0055  Hz  standard  frequency  bandwidth 

AND  30  DIRECTION  BANDWIDTH 
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ENERGY  DENSITY /O 


ALL  DIRECTIONS 


FEBRUARY  1975 


CENTRAL  PERIOD  (SECONDS) 


Table  8:  SPECTRAL  ENTRIES  FOR  A 50-KN0T  FULLY  ARISEN  SEA 

For  a given  location  and  synoptic  time 
0.0055  Hz  and  30  bandwidths 


256 


Figure  1:  CUMULATIVE  SPECTRAL  ELEMENT  STATISTICS 
Curves  are  energy  density  (coded)  per  0.0055  Hz  bandwidth 

AND  30  DIRECTION  BANDWIDTH 
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